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ABSTRACT 

Historically, emission lines have been considered a valuable tool for estimating the 
bolometric thermal luminosity of the accretion flow in AGN, Lboi- We study the re- 
liability of this method by comparing line strengths to the optical/UV continuum 
luminosity of SDSS DR7 radio quiet quasars with 0.4 < z < 0.8. We find formulae for 
Lboi as a function of single line strengths for the broad components of H/3 and Mg II, 
as well as the narrow lines of [O III] and [O II]. We determine the standard errors of 
the formulae that are fitted to the data. Our new estimators are shown to be more ac- 
curate than archival line strength estimations in the literature. It is demonstrated that 
the broad lines are superior estimators of the continuum luminosity (and Ltoi) with 
H(3 being the most reliable. The fidelity of the each of the estimators is determined in 
the context of the SDSS DR7 radio loud quasars as an illustrative application of our 
results. In general, individual researchers can use our results as a tool to help decide 
if a particular line strength provides an adequate estimate of Lhoi for their purposes. 
Finally, it is shown that considering all four line strength, simultaneously, can yield 
information on both Ltoi and the radio jet power. 
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The primary measure of the strength of an active galac- 
tic nucleus (AGN) is its bolometric luminosity, Lboi, the 
broadband thermal luminosity from IR to X-ray. The char- 
acteristic signature of the thermal component is the "big 
blue bu mp" a large blue/UV e xcess in the spectral con- 
tinuum (jSun and MalkanI Il989l l. Typically, one does not 
have complete frequency coverage of the broadband con- 
tinuum and Lboi can only be estimated. The most com- 
mon estimators involve single point rest frame optical or 
UV lumino s ity as an approximation to the big blue bump, 
iKaspi et all (l200Cl|). or t he lu minosity of the optical contin- 
uum as m "i Miller et aJ (l992'). However, there are many in- 
stances in which the optical/UV continuum is not directly 
observable or is contaminated with other sources of flux. 
Thus, astrophysicists need a surrogate for the continuum 
luminosity. The two most prevalent situations that are en- 
countered are flrstly, blazars in which the synchrotron opti- 
cal/UV emission from the relativistic jet is Doppler boosted 
to a brightness that either swamps the quasar thermal emis- 
sion or contributes an unknown fraction of the total ob- 
served optical/UV flux. The other common occurrence is 



when dusty molecular gas (e.g., the "dusty torus") ob- 
scures the thermal emission produced by the quasar accre- 
tion flow from our view. The latter case is believed to be 
representative of NLRGs (narrow line radio galaxies) and 
Seyfert 2 galaxies (|Antonuccilll993l ). Historically, broad line 
strengths have bee n used as a surrogate for continuum lumi- 
nosity in blazars jCao and Jiang' 2001; 'Celotti and Fabian| 
1993i:iCelotti et alill997: GhisoUini ct al 20,11; Gu ct al 20oi; 
Maraschi and Tavecchid l2003l : IWang et all 120041 '). For NL- 
RGs, the broad line region is not visible by deflnition and 
narrow emission line strengths hav e been used as a surro- 
gatc for t he contin uum luminosity l|Rawlings and Saunders! 
, 1991 ; Willott et ajfiggg ). 

These relationships between line strength and Li,oi have 
never been carefully calibrated or scrutinized for their ac- 
curacy. This Letter intends to do both with a carefully 
selected, large sample of SDSS DR7 radio quiet quasars 
(0.4 < z < 0.8) in which two prominent narrow lines ([OUT] 
A5007 and [OH] A3727) and two broad lines H/3 A4861 and 
Mg II A2798 are observable for all objects. This allows us 
to compare the various derived estimators in the context of 
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Figure 1. Scatter plots of the logarithm of the line strengths as a function of the logarithm of Lbol = ISLcont, where Lcont is the 
optical/UV continuum luminosity from 5100 A to 3000 A. The fit to the SDSS DR7 data is given by the black line. The accuracy of the 
linear representation of the data is given by the coefficient of determination, R,^, that is displayed on each plot. An archival estimator 
from the literature is shown as a red line. Notice that the broad component of H/3 is an excellent representation of Lcont at the top 
left. The next best estimator is based on the broad component of Mg II in upper right hand frame. The [OIII] A5007 based estimator is 
clearly superior to [Oil] A3727 based estimator in the bottom row. 



a single sample of objects, which removes the uncertainties 
associated with sample selection biases if a different set of 
objects is used for the calibrations of each individual line 
strength based estimator. In section 4, we consider these es- 
timators in the context of the radio loud subsample of the 
SDSS DR7 quasars (0.4 <z< 0.8). 



1 SAMPLE SELECTION 

In order to determine the dependence of line luminosity on 
the thermal continuum in a quasar, we constructed a sam- 
ple of SDSS DR7 radio quiet quasars with 0.4 < z < 0.8 
and spectra with a median S/N > 7, this yielded 10069 
AGNs. Long slit spectroscopy of radio loud AGN often show 
strong regions of narrow line emission on scales as large as 
~ 100 kpc that t e nd to be aligned with the jet directio n 
l|lnksip et all |2002| : iBest et all [2OO0I : iMcCarthv eral Il995l ). 
The magnitude of this cont ribution to the nar row line lu- 
minosity is largely unknown jWillott et aj[l999l ). Therefore, 
we segregated out the radio loud quasars because there is a 
concern that jet propagation can enhance the line strengths. 
The SDSS DR7 data was cross-referenced to the FIRST 
data base. The radio loudness, R, is usually defined as a 
5 GHz flux density 10 times lar ger than the 4400A flux den- 



sity, R = 'S'sGHz/S'^^jjij^ < 10 (|Kellermann fc Paulinv-TothI 



Il969h . So any FIRST fiux density detection or upper limit 
that implied R < IQ was considered radio quiet with 1.4 
GHz flux density used instead of 5 GHz flux density. We 
also eliminated the low ionization broad absorption line 
quasars from our sample since they are known to have 
anomalously weak [OIII]A5007 and [OII]A3727 emission lines 
jBoroson and GreenI [1993 : IZhang et aill2010l '). In the end, 
there were 6904 radio quiet sources remaining in our sample. 

The spectral data of H/3 and Mg I I regime were sep - 
arately reduced u sing the procedures of iDong et all (|2008l ): 
IWang et all (|2009h which see for further details, and we will 
only briefly outline it here. To measure the H/3 line and the 
[OIII]A5007 line, a single power-law fit to the optical contin- 
uum in the restframe wavelength range 4200 — 5600A was 
obtained taking into account contributions from both broad 
and narrow Fe II multiplets that were modeled using the I 
Zw 1 templates provided by Vcron-C ettv et al.l (2004). The 
H/3 emission lines are modeled as multiple Gaussians: at 
most four broad Gaussians and one narrow Gaussian with 
FWHM < 900 km for H/3, and one or two Gaussians for 
each [OIII]A4959,5007. To measure the Mg II line, we fit the 
continuum from the several continuum windows in the rest- 
frame wavelength range 2200-3500A, after subtracting the 
estimated UV Fell contr i butio ns by the Fe II template gen- 
erated by iTsuzuki et all 1I2OO6I) and the Balmer continuum 
using the method by [Dietrich et U (|2002l ). Each of the two 
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Mg II doublet lines is modeled with one broad five-parameter 
Gauss-Hermite series component and one single Gaussian 
narrow component. Furthermore, the broad components of 
the doublet lines are set to have the same profile. The nar- 
row components are set to FWHM < 900 km and flux < 
10% of the total Mg II flux. Finally, we fit the spectrum in 
the [Oil] regime using one Gaussian for [Oil] emission and 
a single power-law for optical continuum in the restframe 
wavelength range 3600-3800A. This process resulted in five 
pieces of relevant information for every quasar: 

• The optical/UV continuum luminosity from 5100 A to 

3000 A, Lcont 

• The luminosity of the broad component of H/3 A4861, 

• The luminosity of the broad component of Mg II A2798, 

LMgll 

• The luminosity of the narrow line [OIII] A5007, Loiii 

• The luminosity of the narrow line [Oil] A3727, Loii 



2 LINE STRENGTH FITS TO CONTINUUM 
LUMINOSITY 

The blue continuum luminosity is the most basic signature 
of the thermal emission from the quasar, so it is the most 
commonly used quantity for estimating Li,ai- Perhaps the 
most po pular bolometric c orrection is the simple one pro- 
posed bv lKaspi etal HqoO), Lboi « 9ALa(5100A). Clearly, 
using a portion of the optical/UV continuum is more accu- 
rate than a single point and we have that at our disposal. 
Leant ■ For example, a single point could lie in a noisy end of 
the SDSS spectrum. The average spectral index in our sam- 
ple of 6904 quasars, from 5100 A to 3000 A, is = 0.54, 
where av is defined in terms of the spectral lu minosity as 
Li, ^ This spectral slope implies that the lKaspi et all 

l|200(]D bolometric correction can be expressed as 



15Lc 



(1) 



Equation (1) is relatively insensitive to the choice of a, with 
only a few percent change in the constant of proportionality 
as a varies from 0.45 to 0.7. 

Figure 1 shows how well the line strengths represent 
the continuum luminosity with the factor of 15 bolometric 
correction from equation (1). The linear fits in log-log space 
are given by equations (2) - (5) with the standard errors for 
the intercept and the slope 

\og{Ltoi) = 12.32 ± 0.20 + (0.78 ± 0.01) log[LH;3] , (2) 
log(L6„,) = 16.76 ± 0.26 -f (0.68 ± 0.01) log[LM9i/] ,(3) 
log(L6„0 = 26.50 ± 0.32 + (0.46 ± 0.01) log[Lo/ii] , (4) 
\og{Lboi) = 33.96 ± 0.33 + (0.29 ± 0.01) log[Loi/] . (5) 

The coefficients of determination, R'^, are shown on the 
plots. The H/3 fit is the best, followed closely by Mg II 
with the narrow line fits considerably wo rse. Figure 1 con- 
firms the implications of ISimpsonI (|l998l ) that Loiii rep- 
resents Lcont more reliably than Loii does. We also plot 
archival estimators of Ltoi from the literature to show the 
improvement obtained by these more rigorous ca librations 
of the data. It appears that the iKaspi et all (|2000l ) normal- 
ization is considerably higher than what has been used for 



line strength based estimates in the past. A 5100 A estima- 
tor based on the composite of IPunslv and Tingavl l|2006h , 
Lboi ~ 7.4ALa(5100 A), is closer, but still considerably larger 
than the nor malization of t he arc hival line strength based es- 
timates. The Miller et all ll 19921 ') estimator (which is about 
65% of the iKaspi et all (|2000l ) value) seems to match the 
normalization of the line strength based estimates the best. 
However, the poor archival fits are not just a consequence 
of the normalization, but the slop e of the fits as wel l. If the 
normalization were adjusted, the IWang et all (|2004l ) fit for 
H/3 would be fairly accurate. 

There is no unique criteria for the quality of the esti- 
mator and this must be determined by each researcher per 
their requirements. We choose a "factor of two" as a figure 
of merit (i.e. 0.5 < Lostimatcd/^cont < 2, where Lostimatcd 
is the estimated continuum luminosity) for illustrative pur- 
poses. The H/3 based estimate is excellent with 97.4% of 
the estimates accurate to within a factor of two. The other 
estimators are accurate to within a factor of two in the fol- 
lowing order, Mg II, [O III], [O II]; 90.1%, 80.3% and 73.9% 
of the time, respectively. The results are presented graphi- 
cally in Figure 2 as the probability distribution (histogram) 
of the ratio of estimated luminosity to measured contirmum 

luminosity, I/cstimatcd/-f'cont ■ 



3 APPLICATION TO RADIO LOUD QUASARS 

There are 2461 quasars in the SDSS DR7 sample, 0.4 < z < 
0.8, with a FIRST detection and radio loudness, R > 10. We 
use this sample as a test-bed for the application of the esti- 
mators that were derived in the last section. We can check 
our calibrations of the bolometric correction by applying 
Equations (2) to (5) to each line individually and comparing 
the resultant Lestimated to Lcont. The results are presented in 
Figure 3 as the probability distribution (histogram) of the 
ratio of estimated luminosity to measured continuum lumi- 
nosity, I/cstimatod/icont- Our "factor of two" figure of merit 
(i.e. 0.5 < icstimatod/icont < 2) iudlcatcs that the H/3 based 
estimate is still excellent with 94.9% of the estimates ac- 
curate to within a factor of two. The other estimators are 
accurate to within a factor of two in the following order, Mg 
II, [O III], [O II]; 86.2%, 75.1% and 67.5% of the time, re- 
spectively. Figure 3 shows that estimators based on the [Oil] 
and [OIII] emission line strengths have a propensity to over 
estimate the continuum luminosity in radio loud quasars. 
This is in contrast to the radio quiet quasars in Figure 2, 
where the more errant estimates are equally likely to occur 
above or below "1". In other words, [O III] and [O II] are 
often much stronger than expected from the continuum lu- 
minosity alone. The implication is that the excess narrow 
line strength is produced by the aforementioned excited gas 
that is often observed aligned with the radio jets and is 
therefore likely t o be energized by the jet (jlnksip et alll2003 : 
iBest et alllioool : iMcCarthv erail ll995h . We also performed 
the same exercise with a cutoff at 7? > 20. The results were 
very similar to those above, indicating that the definition of 
radio loudness is a negligible factor in this analysis. 

We give an example of how the line strengths can be 
used to analyze radio loud A GN, the blazar 3C 216 . We 
retrieve the line strengths from iLawrence et al.l l|l996l ). For 
H/3, Mg II, [O III] and [O II], from equations (2) - (5), these 
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Figure 2. The histograms represent the distributions of the ratio of the estimated continuum flux to the actual continuum flux, 
^estimated /^cont for each of the emission lines in the context of the SDSS DR7 radio quiet subsample. The smallest dispersion about 1 
(the most accurate estimator) is H/3 in the upper left hand frame. The superior H/3 fit is shown in black in the other three frames for 
the sake of comparison. 



line strengths yield the following estimates for Lboi, 8.12 x 
10"'' ergs/sec, 1.23x10"^ ergs/sec, 7.24x10'''^ ergs/sec, 1.32x 
10''^ ergs/sec, respectively. 3C 216 is argued to be the one 
of the most kinetically dominated known qua sars, with a je t 
power (kinetic luminosity) to Li,oi ratio > 10 l|Punslvll2007l ). 
The broad line estimates are in good agreement, but the 
narrow line estimates are much larger. The results of this 
paper indicate that the broad lines are therefore excellent 
estimators of the continuum luminosity and the narrow lines 
are dominated by a jet contribution. 



4 CONCLUSION 

We have derived line strength based estimators for Lboi for 
radio quiet quasars from the broad components of H/3, Mg 
II, and the narrow lines, [O III] and [O II] in Equations (2) 
-(5). The strength of the broad component of H/3 is a supe- 
rior estimator of Lcont for either radio quiet or radio loud 
quasars since it is accurate to within a factor of two, 97.4% 
and 94.9% of time, respectively. The next best line strength 
based estimators are in order of accuracy, Mg II, [O III] 
and [O II]. We applied our results to radio loud quasars 
and found strong evidence that the narrow line based esti- 
mates are often skewed by what is likely a strong jet induced 
contribution. It was also demonstrated that using all four 
line strengths in tandem can be a useful diagnostic tool for 



studying the jet power, environment and accretion power in 
AGN. 

The results of this Letter are applicable to blazars, but 
perhaps not to NLRGs and Seyfert 2 galaxies. There is com- 
pelling evidence that sometimes the same gas that can at- 
tenuate the nucleus in these sources can als o attenuate the 
narrow line emission (|Mulchaev et aj 1 1994) : iKraemer et all 
I2OIOI ). In general, without additional information, one does 
not know if the narrow line region in a particular Seyfert 2 
galaxy or NLRG is attenuated or not, thus the [O III] and [O 
II] based estimators are not reliable in isolation. We showed 
in section 4 that the narrow line estimators are actually use- 
ful in the blazar context, where attenuation is not an issue. 
In general, narrow line based estimators are most useful if 
complemented by other information such as broad line lumi- 
nosity or the actual continuum luminosity. For NLRGs and 
Seyfert 2 galaxies, the IR continuum luminosity is a valu- 
able complement to the na rrow line strengths since it is not 
believed to be attenuated jFernandes et alll2010l : lOgle et all 
I2OO6I ). 
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Figure 3. The histograms represent the distributions of the ratio of the estimated continuum flux to the actual continuum fiux, 
^estimated /^cont for each of the emission lines in the context of the SDSS DR7 radio loud subsample. The smallest dispersion about 1 
(the most accurate estimator) is produced by the H/3 based estimate in the upper left hand frame. The superior H/3 fit is shown in black 
in the other three frames for the sake of comparison. 
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